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Abstract 

For  a  system  of  a  multilevel  atom  interacting  with  a  single-mode  cavity 
field,  squeezing  of  the  field  is  studied  numerically.  It  is  found  that  the 
squeezing  effect  becomes  stronger  when  the  atomic  level  number  increases,  but 


tends  to  saturate  when  the  number  exceeds  twenty. 
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I.  Introduction 

The  interaction  of  an  atom  with  cavity  fields  has  stimulated  great 
interest  because  many  interesting  and  novel  phenomena  have  been  discovered. 
In  addition  to  quantum  collapse  and  revival  (Eberly  et  al  1980),  the  atomic - 
field  system  has  been  shown  to  exhibit  a  subPolssonlan  photon  distribution 
(Puri  and  Agarwa  1986),  squeezing  (Meystre  and  Zubairy  1982)  and  stationary 
atomic  inversion  (Savage  1988).  Recently,  Li  and  Zhu  (1985)  and  Kozlerowski 
(1986)  studied  the  dynamics  of  an  M-level  (M-l)-mode  system  with  the  atom 
having  a  common  upper  level.  Most  recently,  the  present  authors  (Li  et  al 
1989)  examined  the  effects  of  the  atomic  level  number  on  quantum  collapse  and 
revival  and  photon  statistical  properties  in  the  interaction  of  a  multilevel 
cascade  atom  with  single>mode  cavity  field.  In  the  present  paper  we  study 
squeezing  effects  for  an  M-level  atom  interacting  with  a  single-mode  cavity 
field.  First  we  present  the  formalism  for  the  problem.  Then  we  proceed  to 
study  fluctuations  of  the  field  mode  and  examine  effects  of  the  degenerate 
multiphoton  processes. 

II.  Formulation 

We  consider  a  cascade  atom  with  M  energy  levels  as  shown  in  Fig.  1.  The 
total  Hamiltonian  is,  in  the  rotating-wave  approximation, 

M  M-1 

H  -  >(Oa^a  +  ^  ^  ^i^i+l^i*  +  h.c.  ,  (1) 

i-1  i-1 

where  the  operator  a^  creates  a  photon,  a|  creates  an  atom  in  the  i-th  level, 
and  are  the  atom- field  coupling  constants.  Our  discussion  will  be  limited, 
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for  simplicity,  to  processes  involving  only  a  single -photon  transition  between 
adjacent  atomic  levels. 

We  define  the  total  excitation  number  operator  as  (Li  et  al  1989) 

M 

N  —  a^a  +  ^  ^^i^i  ’ 

i-1 

which  is  the  sum  of  the  photon  number  and  atomic  excitation  number  operators. 

A  A 

It  can  readily  be  verified  that  [H,N]  -  0,  i.e.  ,  N  is  a  constant  of  motion. 
If  we  are  only  interested  in  the  particular  case  in  which  the  atom  is 
initially  in  its  higher  level,  then  we  need  only  consider  that  part  of  the 
Hilbert  space  involving  state  vectors  corresponding  to  N  >  M,  where  N  is  the 

A 

eigenvalue  of  the  operator  N.  This  is  a  consequence  of  the  fact  that  N  is  a 
conserved  quantity. 

Consider  the  state 


I  i ,  n>  -  I  i>  I  n> 


(3) 


in  which  there  are  n  photons  in  the  field  and  the  atom  is  in  the  i-th  state. 
The  subspace  corresponding  to  N  —  M  +  n  is  spanned  by  the  state  vectors  |M,n>, 

|M-l,n+l> .  |i,n+M-i> .  |l,n+M-l>,  where  n  is  the  photon  number  in  the 

field  and  can  take  any  positive  integer  value.  Since  state  vectors  in 
different  subspaces  are  necessarily  orthogonal,  the  total  Hamiltonian  can  be 
diagonalized  in  every  subspace. 

An  arbitrary  state  in  the  subspace  corresponding  to  N  can  be  expressed 


as 


Availability  Codes 


□  □ 
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'V-  2  ■ 


where  the  expansion  coefficients  C.  „  .  satisfy  the  s tationary- state 


Schrodinger  equation 


^  ^^ii^ii'  “ii'  ‘  ^^ii' ^^i' ,n+M-i' 


The  matrix  elements  in  (5)  are  given  by 


H..  -  <i,n+M-i|H|i,n+M-i>  -  (n+M-l))tf)  + 


A.  -  a 


^ii'  “  <i-"+M-i|H|i' ,n+M-i'> 


-  >j^,yn+M-i'  +  Aj^,_j^yn+M-r+l  ^  "  i'  -(^b) 


Here  we  have  defined  the  detuning  par£uneters 


-  Htl  -  .  i  -  1.  2 . M-1  .  (6) 


The  energy  eigenvector  «nd  the  corresponding  eigenvalue  are 


determined  by 


hU  >  -  E  U  > 
'  no  no'  no 


(7a) 
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M 


\<f>  > 

'  na 


i,n+M-l 


i,n+M-i> 


i-1 


(7b) 


The  subscript  a  in  (7a)  labels  the  eigenstates  in  the  subspace  in  question. 

If  the  atom  is  assumed  to  be  initially  in  the  M-th  level,  then  for  all 
time  t  >  0  the  system  can  only  be  found  in  subspaces  with  n  >  M,  There  are  M 
eigenstates  of  the  type  (7b)  in  every  subspace.  Hence  the  orthomormality 
condition  is 


M 

i-1 


i,n+M-i  i,n+M-i 


a'  a 


(8) 


and  the  completeness  relation  is 

«  M 

n-0  <7-1 


We  can  now  solve  the  equation  of  motion  for  the  density  matrix.  The 
reason  why  we  use  the  density  matrix  Instead  of  the  state  vector  formulation 
is  that  it  is  easier  for  this  to  be  generalized  to  treat  the  complicated 
problem  of  a  mixed  atomic  initial  state,  which  we  shall  study  next.  Besides, 
once  the  Hamiltonian  H  has  been  diagonalized,  the  elements  of  the  density 
matrix  can  be  obtained  in  a  straightforward  fashion  without  any  increase  in 
complexity.  In  fact,  even  if  one  uses  the  state  vectors,  one  still  has  to 
diagonalize  H.  A  density  matrix  element  between  energy  eigenstates  satisfies 
the  equation 
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dt  ^na.n'a' 


M  ^n'<7' ^^na.n'a' 


(10) 


whose  solution  can  be  put  in  the  form 


^na.n'a'^^^  "  ^no.n'a' ® 


-~  (E  -E  ,  ,)t 
)<  no  n  a 


(11) 


The  initial  matrix  element  in  (11)  is  given  by 


(12) 


where  C  ^  is  the  inverse  matrix  of  the  transformation  defined  in  (4)  and  p  , 

nn ' 

is  the  initial  density  matrix  element  of  the  field.  The  last  step  is  from  the 
inverse  transformation  of  (7b),  namely, 


M 


n,M>  - 


(C*^)^  \<k  > 

^  '^<7,n+M-a'’^na 


(13) 


(7-1 


If  the  Initial  density  matrix  is  known,  the  mean  values  of  all  dynamical 
variables  of  the  atom* field  system  can  be  calculated.  The  mean  occupation 
probability  of  the  j-th  level  is 


Pj  -  Tr(pA 


■jS>  -  \  \ 


P  ,  0^^  0  0  I  ^1^4  1^  ^ 
nV' *  j  y^na 


na  n'a 


-} 


n,a,a' 


(14) 
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It  is  easily  verified  that 


^na.nor'  /  j,n+M-j  j,n+M-j 


n<7 ,  n<7 


n.CT.CT' 


where  we  have  made  use  of  (8)  and  (9). 


As  usual,  we  introduce  the  slowly-varying  complex  amplitudes  for  the 


field, 


(16a) 


d2  -  -  ate-'<«"*^>l 


(16b) 


These  operators  satisfy  the  commutation  relation 


[d^.d2]  -  J  . 


For  simplicity,  in  this  paper  we  only  consider  the  case  of  ^  -  0,  for  which 
maximum  squeezing  is  expected,  at  least  for  the  case  where  the  atom  starts 
in  a  definite  state  (Meystre  and  Zubairy  1982).  The  corresponding  expressions 
for  the  fluctuations  are  then 


(Adj^)^  -  J[2<a^a>  +  1  +  <a^>e^^®*^  +  <a^^> 


^  iOt  ^  t  -int.,^2, 
-  <ae  +  a'e  > 


(18a) 


(Adj)2  .  l(2<.t.>  t  1  -  <.'>.2“'  -  <.t^> 
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^  ^  int  t 

+  <ae  -  a ' e  >  ] 


(18b) 


where 


<a  a>  -  Tr<pa 


p  ,,<^,,a'a^  > 

Tuj,n'a'  n'a'  '  ’  tw 


na  n'a 


-  )  (n+M-i)  )  p  ,  (0)  C°  ^  .  C°'  „  .  cos{E  -E  ,)t)  , 

l_  ^  L  n^.n  CT  j,n+M-j  j  ,n+M-j  na  na' 


<ae  +  a'e  >  -  Tr[p(ae  +  a'e  )] 


-  2  )  )  Vn+M-i+l  )  ,  ,  (f'  ,  p  /  ,(0) 

i  l_  n+M-i, i  n+M-i+l,i  '^na,  (n+l)a' ' 


X  cos[ ( 


E  -  E. 
na  ( 


+  n)t]  , 


(20a) 


^  iOt  t  -iOt  _  ,  .  iOt  t 

<ae  -  a'e  >  -  Tr(p(ae  -  a'e  )] 


E  -  E, 


P  /  \  »  (0)  sl"l  (■ 

na,(n+l)a'' 


+  n)t]  (20b) 


^  2  2int  ^  t2  -2int  -  r  /  2  2int  ^  12  -2iOt, 

<a  e  +  a'  e  >  -  Tr[p{a  e  +  a'  e  ) 


f(n+M-i+l)(xH41-i+2) 


'n+M-i,  t 
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E  -  E 

.  .  ng  ('n4-2)g'  , 

^  •  P  /  ns  I  (0)  cos[(  w  +  2n)t] 

n+M- 1+2,1  n<T,(n+2)CT'  )< 


(20c) 


They  satisfy  the  uncertainty  relation 

Ad-  ‘Ad.  >  '7  .  (21) 

1  2  ~  4 

When  either  of  them  satisfy  the  condition 

(Ad.)^  <\  .  i  -  1,  2  (22) 

we  say  that  the  field  is  in  a  squeezed  state. 

For  a  given  n,  Eq.  (5)  can  be  solved  for  the  energy  eigenvalues  and 
eigenfunctions.  With  the  initial  density  matrix  elements  detexrmined  by  (12), 
we  can  calculate  the  field  fluctuations  and  occupation  probabilities  of  any 
atomic  level  at  any  time  t  >  0. 


III.  Results  and  discussion 

For  simplicity  we  assume  that  ”  ^2  ”  '  "  '  "  ^M-1  ”  ^  that  the 
resonant  interaction  condition  is  satisfied  (Li  and  Zhu  1985,  Kozierowski 
1986,  Li  et  al  1989),  i.e.,  Aj^  -  0.  In  our  numerical  calculations  we  take  )(  - 
1  and  employ  1/A  as  the  unit  of  time.  For  the  initial  conditions,  we  consider 
the  atom  to  start  from  the  uppermost  level  |M>  and  the  field  mode  to  be  in  a 
coherent  state.  Thus  we  have 


n'  .n  - 
^nn'  "  yn!n' ! 


(23) 
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where  n  - 


is  the  initial  mean  photon  number. 


The  relation  between  the  field  squeezing  and  the  initial  field  intensity 


n  is  nonlinear.  For  a  given  energy  level  number  of  a  cascade  atom,  there 

exists  a  value  of  n  corresponding  to  the  strongest  squeezing  or  the  smallest 

2  2 
(Ad^)  (Zhu  et  al  1988),  which  hereafter  is  referred  to  as  ( 

maximum  squeezing  corresponding  to  different  level  number  M  for  some  cases  can 


be  seen  from  Fig.  2.  Some  of  the  results  are  summarized  in  Table  1,  from 
which  it  is  observed  that  the  field  squeezing  increases  with  the  increase  of 
the  cascade  level  number.  It  is  also  noted  that  the  larger  the  value  of  M, 


the  later  the  field  enters  squeezed  states  and  the  longer  it  remains  squeezed. 


However,  on  the  other  hand,  we  find  that  there  is  saturation.  For 

2 

example,  when  M  changes  from  2  to  3,  decreases  from  0.2023  to  0,1730; 

2 

when  M  changes  from  8  to  10,  decreases  from  0.1309  to  0.1243;  when  M 

changes  from  21  to  23,  decreases  from  0.1090  to  0.1077.  Detailed 

2 

behavior  is  shown  in  Fig.  3.  In  brief,  for  larger  M,  much 

less.  The  variance  tends  to  saturate  to  some  non-zero  value. 


Figure  4  gives  an  example  of  the  time  evolution  of  the  occupation 
probability  for  the  atom  in  the  uppermost  state.  Comparing  with  Fig.  2,  we 
can  see  that  whenever  the  field  enters  the  squeezed  state  from  its  initial 


coherent  state,  the  probability  of  the  atom  lying  in  the  uppermost  state  is 
close  to  its  minimum.  The  opposite  situation  occurs  when  the  field  leaves  the 
squeezed  state.  Thus  we  can  conclude  that  squeezing  effects  have  something  to 
do  with  the  absorption  process,  which  is  in  agreement  with  the  results 
reported  by  Shumovsky  et  al  (1987). 

Ue  have  discussed  above  the  effects  of  the  quantities  M  and  n  or  z  on 
the  squeezing.  As  for  the  effects  of  the  other  quantities  and  we  find 
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that  increasing  detunings  would  weaken  the  squeezing,  and  there  would  be 
stronger  squeezing  when  the  coupling  constants  X’ s  are  equal. 
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Table  1:  Maximxun  squeezing  under  different  conditions 


M 

n 

maximum  squeezing 

2 

6 

19% 

4 

7 

36% 

6 

9 

43% 

10 

11 

50% 

20 

15 

56% 

14 


Figure  Captions 

1.  Schematic  energy  level  diagram  of  the  M-level  atom. 

2  2 

2.  Evolution  of  the  variances  (Ad^)  (curve  1)  and  (Ad2)  (curve  2) 

-  2,  ii  -  6;  (b)  M  -  8,  n  -  10;  (c)  M  -  19,  n  -  15. 

2 

3.  The  minimum  variance  (Ad,)  .  versus  level  number  M. 

1  min 

4.  Evolution  of  the  highest  level  occupation  probability  for  M  - 


(a)  M 


8  and  n  - 
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